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ABSTRACT 



Ths PERT/COST system, now in current use tequlr- $ several simplify 
ing assumptions in order to permit the use of an&iytic methods for deter- 
mining predictions of critical path activities and project duration. As 
a result of these assumptions, the predictions may involve significant 
error, and will always be optimistic* The Monte Carlo technique of net- 
work calculation does not require these assumptions and* hence 9 is 
capable of yielding more accurate predictions and providing more useful 
information. This technique is discussed in detail ■ A more flexible 
probability density model for activity times is introduced* and a re- 
source allocation technique based on the probability that an activity 
will be on the critical path is developed* Finally* an application of 
PERT network theory to military operational planning is described* 
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NOTATION 



Symbol 


Introduced 
on page 


Meaning 


* 


9 


Indicates a random variable 


t* 


9 


Time required to complete an activity not yet 
completed 


t 


7 


A specific value of the random variable t* 




11 


Expected value, or mean, of t* 


t s 


25 


Time allowed by schedule for completion of an 
activity 


T* 


12 


Elapsed time till the occurrence of an event, or 
until completion of an activity 


t e 


12 


Expected value, or mean, of T* 


T S 


26 


Elapsed time at which an event is scheduled to 
occur 


Tl * 


13 


Latest elapsed time an event can occur and not 
delay the scheduled completion of the project, or 
the latest time an activity can be completed and 
not delay the project 


Tl 


13 


Expected value, or mean, of Tj* 


Se 


25 


Earliest elapsed time an activity/event can be 
scheduled for completion 


S L 


25 


Latest elapsed time at which an activity/event can 
be scheduled without causing a change in t s for 
some subsequent activity 


K* 


13 


Slack for an activity /event , = T^* - T* 


k e 


14 


Expected value, or mean, of K* 


</( *) 


10 


Standard deviation of the random variable in 
parentheses 


d' 2 { *) 


12 


Variance of the random variable in parentheses 


D* 


14 


Duration of the entire project 
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NOTATION (cont’d) 



Introduced 
Symbol on page 

D„ 8 
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14 



Meaning 

Elapsed time to the Directed Completion Date for 
the project 

Probability an activity/event will be on the 
critical path. The Criticality Index. 

Probability an activity/event will be completed 
on schedule 
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Probability that the project will be completed 
on or before the DCD 



DCD 

a 
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10 



Directed Completion Date for the entire project 

Optimistic estimate for t*, the smallest amount 
of time the activity can reasonably be expected 







to consume 


M 


10 


Most likely estimate for t*. The elapsed time at 
which the activity is most likely to occur. 


P 


10 


Pessimistic estimate for t*. The greatest amount 
of time the activity can be expected to consume. 


X* 


41 


The standard Beta variable 


m 


41 


The mode of x* 
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Introduction. 



The PERT system was devised by a study group under the direction of 
the Special Projects Office, Bureau of Naval Weapons, Navy Department, 
in order to develop a methodology for providing the management of the 
Fleet Ballistic Missile (Polaris) program with an information reduction 
system for program monitoring and evaluation. By this system, management 
was to be continuously apprised of progress to date for the program as a 
whole, and be furnished valid predictions as to outlook toward accom- 
plishing program objectives, \l\ 

The system, as developed, was a giant stride in management tech- 
nology, and is credited with a major contribution toward the rather 
phenomenal success of the FBM development program* Not only did the 
PERT system accomplish the objectives already mentioned, but it allowed 
management to predict those activities in the development project whose 
completion time would have a direct effect on the overall project dura- 
tion. These critical activities could then receive appropriate mana- 
gerial attention. 

The PERT system has since been employed extensively in industry, 
and has been made a standard procedure for monitoring research and 
development projects under the cognizance of the Federal government. 
Industry acceptance of the system has been, in general, enthusiastic and 
widespread. 

The system has been recently extended to include a cost control 
feature, PERT/COST, which has now been designated the standardized sys- 
tem for monitoring federally sponsored R&D projects. 

Another system, the Critical Path Method, has paralleled PERT in 
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development. It is similar in many respects to PERT, lacking PERT's 
stochastic representation of time, but incorporating a correlation be- 
tween time and cost, enabling management to schedule optimally. The 
Department of Defense and NASA have incorporated the desirable features 
of both systems into PERT/COST. 

The success of PERT and Pert-type systems has stimulated the inter- 
est of professionals in the fields of mathematics, management science, 
and operations research. The professional journals have featured many 
papers on the relative merits of the sys terns 9 proposals for refinement 
and extension of the systems, consolidation and integration of the sys- 
tems, and the like. It is the purpose of this paper to consolidate many 
of these ideas, and to propose some refinements and extensions to the 
system, designed to extend its usefulness and increase the accuracy of 
predictions about the projects, and to generate specific recommendations 
for control of time and cost. 

We shall first briefly describe the basic PERT system, critically 
discuss certain features of the system, and tender proposals by which 
the basic system may be improved. In particular, we shall discuss the 
Monte Carlo method of calculation of the PERT network, and show how this 
method overcomes many of the deficiencies inherent in the analytic method. 

Next we shall describe the standard PERT/COST system currently being 
implemented by DOD and NASA. We shall follow this with our proposal for 
an integrated project control, scheduling, and resource allocation tech- 
nique . 

Finally, we shall discuss an extension of the PERT network approach 
to the problem of planning projects involving major future decisions, 
such as military campaigns. 
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2 . 



Conclusions and Recommendations 



As a result of this study, it is concluded that: 

(a) The assumptions employed in the currently standard analytic 
PERT network calculations may induce significant error in the results . 

(h) The standard analytic technique for network calculation fails 
to provide a sufficiently definitive measure for activity criticality. 

(c) A more flexible model for activity completion time, providing a 
range of variance to allow for varying uncertainty in time estimates, is 
needed. 

(d) The Monte Carlo technique does not depend on the assumptions 
which induce the inaccuracy in the analytic technique, hence is capable 
of yielding results limited in accuracy only by the validity of the 
estimates for activity time. In addition, this method yields the proba- 
bility that a given activity will be on the critical path, and other 
useful information not obtained by the analytic technique. The Monte 
Carlo technique can easily provide the flexibility discussed in (c), 
above , 

(e) The Monte Carlo technique requires increased computer running 
time, but the increase is not prohibitive, and is well justified by the 
improvement in quality of the results obtained* 

(f) A resource allocation technique, based on the probability that 
an activity is on the critical path, is feasible, and provides a reason- 
able method of optimum resource allocation. 

(g) The network representation of the complex interrelationships 
between activities in a project is a technique which may well be applied 
to the task of planning military operations. 
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It is recommended that the Monte Carlo technique for PERT network 
calculations be employed in lieu of the analytic technique. A method of 
utilizing this technique is discussed in detail in Appendix I. 
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3. The Basic PERT System 

The basic PERT system is a means of defining the relationship be- 
tween the various activities comprising a project*, and for estimating 
the time required to complete each activity and the project as a whole. 
Use of the system enables realistic schedules to be imposed on the 
project* and assists in controlling the execution of the project. 

By a P r °J ec t we shall mean a collection of tasks, or activities * 
each a necessary step toward the achievement of some final objective,. 
Subsets, called paths , of these activities, are dependent sequentially , 
that is, any activity in the path, other than the beginning activity , 
may not be commenced until its predecessor has been completed. Each 
activity in the path, except the beginning activity*, has a predecessor 
in the path. An example of a project in this sense may be as simple as 
the construction of a house, where the activities are tasks such as lay- 
ing foundation, erecting subfloor, installing rough plumbing, erecting 
framing, installing rough wiring, etc. Or a project may be as complex 
as development of the Fleet Ballistic Missile system, with many hundreds 
of activities involved. 

A network of directed arcs is a very convenient way of representing 
a project. The activities are represented fey directed line segments, 
called arcs, terminating at nodes, called events. (See Figure 1.) The 
time required to complete an activity is associated with the length of 
the arc representing it. Events represent points in time. Several 
activities may terminate at one event. If so, the event is said to occur 
when the last of these terminating activities is completed. At that 
point in time, other activities, whose commencement is contingent upon 
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Figure 1 

A Network Representation of a Project 

the event, may commence. These activities are represented by arcs di- 
rected outward from the event* In the house construction example, the 
installation of the interior wallboard may commence only after the rough 
plumbing and rough wiring are completed* 




Figure 2 

Two Independent Activities Incident to an Event 

A project can always be represented by a network with one beginning 
event and one end event. Fictitious, or dummy 9 activities may be drawn 
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in the network in order to achieve this format,. Such dummy activities 
may consume zero time, or may represent a predetermined amount of sched- 
uled ’’dead” time* It is common practice for the network to be drawn 
with beginning event on the left, and with the end event, representing 
the achievement, of the final objective, at the extreme right of the net- 
work. 

A path may be defined as an unbroken chaiin. of events, with origin 
at the beginning event and terminus at the end event <> The network con- 
sists of many such paths, overlapping, paralleling, and crossing each 

other. 

Activities generally consume time. The exact amount of time, t, 
which will be consumed in completing any one activity is not generally 
known in advance. If one looked at the history of a completed project, 
one could discover the exact amount of time, t } required to complete any 
given activity, and could label each event with its actual time of occur- 
rence. The longest path, with respect to time, from the beginning event 
to the end event could then be found. This path is called the critical 
path , The activities on this path are called critical activities . The 
length of the critical path, the sum, ^t, of the completion times for 
all activities on the path, is the total elapsed time from the beginning 
to the end of the project, and is called the duration , D* It is apparent 
that a small change in t for one of the activities on the critical path 
causes the same amount of change in D, but a small change in t for an 
activity not on the critical path would cause no change in D. Hence it 
is important that project management be able to predict the critical 
path in order to predict D* and rm order to allocate resources 9 and fix 
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schedules, most effectively. 



Any path, other than the critical path, is shorter than the critical, 
path. The difference in lengths of the critical path and any other path 
is called the slack , K, for the n©n~critical path. An activity or event 
may be located on several paths simultaneously. The smallest slack 
associated with any path which includes a given activity or event is the 
value of slack associated with that activity/event. The activities/ 
events on the critical path have zero slack. Any single activity may be 
delayed an amount less than or equal to its slack without affecting the 
duration of the project. 

The term, slack , may have a slightly different meaning if some au- 
thority has imposed a Directed Completion Date , BCD* on the end event of 
the project. If so, we define the Scheduled Duration , D g , as the elapsed 
time from the beginning event to the DCD» Then slack for the critical 
path is (D s - D) . Slack for any other path is (D s * path length). With 
this definition, slack for any activity may be either positive, negative, 
or zero, but all activities on the critical path will have the same 
slack, which will be less than that for any other activity in the net- 
work. 

The PERT system is a method by which predictions of the critical 
path, and slack for any activity* may be made with sufficient accuracy 
to enable management to operate more effectively. 

The Stochastic Ne twork Model . In the planning stage* before a proj- 
ect has commenced, the actual time* t* which will be required to complete 
any given activity, is an unknown quantity. It is unlikely that this 
time could be predicted exactly. We shall call such uncertain quantities 
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random variables , and symbolize them with an asterisk f as t** The corre- 
sponding unstarred symbol will represent a specific value of the variable 
It is possible to set limits, within which t* is almost certain to fall. 
It is also possible to estimate the shape of the probability density 
function governing the random variable 9 t*. The probability density 
function for t* is a curve whose ordinate in a region is a measure of 
the likelihood that t* will fall in that region* To be more precise, 
the area under the curve in a region is the probability that t* will 
occur in that region, in Figure 3 the ratio of the shaded area to the 
total area under the curve is the probability that t* will fall between 
and t2* 

m ^ 



/ 






\ 



X 






-1 

Pigure 3 



A Probability Density Function 



The actual probability density function for t* is difficult, if not 
impossible to obtain, and will be unique for any given activity performed 
in a given time period. Theoretically, if the activity were performed 
repeatedly, under the same conditions, with no learning taking place. 
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then a histogram showing the number of completions occurring in a time, 
t i At, plotted against t in discrete segments , would approximate the 
probability density function. Since this procedure is not possible* we 
must be satisfied with the best estimate we can formulate* using avail- 
able knowledge of the nature of the activity* 

In order to estimate the probability density function of t* for a 
given activity,, the most qualified supervisor in charge of the activity 
makes three estimates of the completion time for the activity. One 
estimate, an optimistic estimate, is the smallest amount of time the 
activity may be reasonably expected to consume* Another, M* is the 
amount of time the activity is most likely to require* Finally, he 
makes a pessimistic estimate, P, which is the greatest amount of time 
the activity can be expected to consume* barring completely unforeseen 
circumstances* P and 0 form the upper and lower bounds* respectively* 
of the variation of t*, and M is the mode of the density function* The 
probability density function is assumed to have the value zero at P and 
0, and to reach a single peak at M* 

Any of several well known mathematical functions exhibit the prop- 
erties we have just assumed for the distribution of t*. The Beta 
function, 

f ( C) ; kCt-a)* (b-t/ (1) 

with parameters a t b,o< , and / , was chosen as the model for the PERI sys- 
tem activity time because it could be specified to fit the constraints 
of the estimates. The three estimates 0, H* and P* specified three of 
the parameters of the Beta function. The remaining parameter was sped 
fied by arbitrarily requiring the standard deviation* efft*}, to equal 
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ILzA . One-sixth of the range is a frequently used estimator for the 
6 

standard deviation of unimodal frequency distributions , hence was con- 
sidered to be a reasonable assumption* For example, the standard normal 
dis trihucion, truncated aft 2.66, has its standard deviation equal to 
1/6 the range. [_6] As a result of the above requirements, one of the 
parameters, say ot , becomes a function of the relative position of M in 
the range between # and P„ Determination of this parameter requires the 
solution of a cubic equation, (see Appendix I). The exact determination 
of the mean, or expected value of t*, symbolized by t e , is a tedious 
process, since t e is a function of a,b 9 oCand /. The following reasonably 
close linear approximation for the mean was adopted as a standards 

E [t*J 5 t £ , t ±M±J_ ( 2 ) 

6 

The probability density function thus assumed has a shape similar 
to that shown in Figure 4, This distribution will be referred to subse- 
quently as the PERT Beta distribution. 




A Typical PERT Beta Density Function 

A raore detailed discussion of the PERT Beta distribution is included in 
Appendix L 
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For purposes of predicting the critical path, the three estimates 
are combined by means of the weighted average formula, (2), for the mean* 
The expected value, or mean, t e , is used as a deterministic value of t* 
in the prediction of the critical path* t e lies between H and the mid- 
point of (#,P), and is separated from M by 1/3 the distance between M 
and the midpoint. In this manner, account is taken of the usual tendency 
for estimators to be optimistic about placing M. The fundamental reason , 
however, for using t e instead of M as an estimator for t* lies in proba- 
bility theory, rather than an effort to correct for this supposed bias* 

In order to determine Slack , two new random variables, T* and Tj*, 
are defined. T* is associated with a given event, and is the elapsed 
time between the beginning event and the occurrence of the given event* 

T* is the sura of the random variables, t** for the activities on the 
longest path to the event. At this point in the development of the 
basic PERT system an assumption is made that the longest path to an 
e vent will be the path whose mean length is greatest . The mean length of 
a path is the sum of the means, t e , of the activities on the path. Then 9 
by the above assumption, the mean, Tg, of T*, is the largest of the mean 
lengths of all paths to the event. Operating on the basis of the same 
assumption, the further assumption is made that T* is approximately nor- 
mally distributed. This assumption is justified since we assume we know 
the longest path to the event, and since this path is the sum of n random 
variables, the Central Limit Theorem of probability theory asserts its 
tendency toward normality as n increases. Due to the assumed independ- 
ence of the activities along this path, the variance, d ^(T^) , of is 
the sum of the variances, of the activities on the path to the 
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event. The distribution of X* is now completely specified by Tp and 
O' 2 ( T* ) in view of the normality of T*. 

We may associate a random variable, X*, with an activity as well 
an event. We define T* for an activity as the elapsed time until the 
activity is completed . Note that T* for an activity is not necessarily 
equal to T* for the event following the activity. T* for that event is 
equal to the maximum of the T*"s for the activities incident to it, that 
is, terminating at the event. 

The random variable, for an event is defined as the latest 

time the event can occur and not delay the scheduled completion of the 

project. Tj* is calculated by subtracting from D s the suns of the t*’s 

for all the activities on the longest path from the event to the end 

event. Again the assumption is made that this longest path will be the 

path whose mean length, longest. Then the mean value of 

written T^> - (D s - ^t e ), where the sum is the maximum of all such sums 

on paths between the given event and the end event. The variance of 
2 9 

Ti* t d (Tj*), = d (t*), along the path with greatest mean length, 

A landom variable, Ty*, may also be associated with an activity, 

Tj* for an activity is defined as the latest elapsed time since the be- 
ginning event at which the activity must be completed in order not to 
delay the scheduled completion of the project. Notice that Tj* for any 
activity is the same as the Ty* for the event following the activity. 

It is apparent that the sum of the two paths used to determine T* 
and iy* for a given event/activity is the longest path from the begin* 
ning to the end which contains that event/activity. From our definition 
of slack, we obtain the following relationship for slack, K*’ 
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( 3 ) 



K* : D s - |t* (D s - Tj*)J ; Tj* - X* 

Then the mean value. Kg » of 5 Jack is the difference, T r ~ . Notice 

that this expression for Kg depends upon the assumptions made in deter- 
mining and Tg« 

The Probability of Meeting a Scheduled Completion Date. We have 
developed the random variable, T*, whose distribution we have specified „ 
by means of simplifying assumptions, to be normally distributed, with 
mean Tg and variance, O' ^(T*) . The probability, Fg, that an event will 
occur before some scheduled completion date is simply where 

T s is the elapsed time from the beginning event to the scheduled comple- 
tion date. In particular, the probability, that the project will be 

completed before the directed date is P [p* < D S "J * Where the duration, 

D*„ is a random variable equal to T* for the end event. These probabili- 
ties can be obtained from a table of values for the normal random 
variable. 

The Predicted Critical P ath is that path whose mean value of slack. 
Kg, is smallest. All events or activities on this path will have the 
same value for Kg, which will be less than the Kg for all other events/ 
activities. 

This completes the description of the basic PERT system which is 
currently in general use. In practice, the network representing the 
project is drawn, events and activities are assigned identification num- 
bers, and are assigned to specific organizations for responsibility. 

The supervisors responsible for completion of each activity then submit 
the three estimates, 0, M, and P* Several computer programs have been 
written which take the activity number, preceding and succeeding event 
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numbers, and the three time estimates as inputs* then perfom the net- 



predetermined ordering* giving Tg 8 T^* Kg* the probability* of meet- 



critical path, and several other paths in increasing order of Kg, Man- 
agement may take appropriate action by reallocation of resources if the 
probability of meeting the overall project BCD is unacceptably small. 

In order to "buy time 11 additional resources are allocated to* and addi- 
tional managerial attention devoted to those activities on the predicted 
critical path, and to those whose mean slack is relatively small. Re- 
sources may be shifted from paths with greater Kg to paths with lesser 
Kj ? t In some cases reorganization of the network is required* placing 
more activities in parallel rather than in series paths* or eliminating 
activities whose products are not entirely essential to the project,, 

The results of original PERT computations may be the basis for 
establishing an overall project schedule. This schedule specifies sched- 
uled completion times* t s , for activities* and scheduled times ©f occur- 
rence* for certain events. 

While the project is in progress* revised inputs are submitted 
periodically, showing percent of completion* and revised estimates. In 
the case of a partially completed activity* these estimates should be 
considerably more reliable than those submitted before work was commenced. 
These inputs are fed to the computer* and another network calculation is 
run, as before, except that completed activities are noted and the actual 
completion time is used in the network rather than an estimated time. 

The results of these runs may necessitate schedule changes* and/or any 



work calculations. 




The output lists the activities in any 



. 2 

ing a DCD„ d *(T*) , t e , and (f“{ t*). The report may list the predicted 



of the actions previously Mentioned. 
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Critique of the Basic PERT System 

The system which we have just described tends to identify and call 
attention to a single path, the Predicted Critical Path, This path, 
though actually uncertain, is deterministically established by virtue c£ 
several assumptions * In most networks, the activities on this "critical 
path' 1 will have probabilities of actually being critical much less than 
unity. Furthermore, other activities, not mi this path, my have proba- 
bilities of being critical larger than some activities on the "critical 
path". Ij&J Hence, the system, by focussing attention on the single path, 
may de emphasize activities which are actually more important. 

The assumptions which were made in determining Tg and lead to 
optimistic results in these figures, and in Kg, By reason of the same 
assumptions, the probability that an event is completed before a DCD is 
always optimistic, 

is the random variable representing the elapsed time to occur- 
rence of an event* Suppose there are n different paths to the event. 

Let the random variables, p j*, p^/*, be the lengths of these n 

paths. Mow T* is the maximum of these n variables, but the maximum of a 
set of random variables is not normally distributed, even if the indi- 
vidual variables themselves are normal* In particular, T* does not have 
the same distribution as the p* with the largest mean, nor does T* have 
the same mean as this p*. Furthermore, the p^* are neither independent 
nor identically distributed, hence the calculation of the distribution 
of the maximum of this set is so complicated as to render it infeasible 
by analytic methods. 

In general, the probability that an event occurs before some 
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specific scheduled date is the probability that none ©f the possible 
paths overrun the scheduled date. By the PERT assumptions, we consid- 
ered only one path, namely that with the greatest mean length. Ob- 
viously, if there are several parallel paths with nearly the same mean 
length, then by considering only the probability that one of the paths 
does not overrun, we are being quite optimistic in our result. 

For example, suppose there are three parallel and independent paths, 
each approximately equally likely to be the critical path* and each hav- 
ing a probability of being less than D s of about Q.5Q* Then P , the 

probability that the project is completed prior to the DCD, is the pr©b~ 

3 

ability that none of the path lengths exceed D s « Hence* a (.5) ® .125. 
The PERT procedure would consider only the path with the largest mean 
length, and would calculate P D s 0.50, For networks with several paral- 
lel paths of comparable length, the error would be quite large , as in 
the example. In networks with only one predominantly long path, the 
error would be negligible. 

In order to obtain a more general prediction of criticality than 
the single path prediction, we need to calculate the probability that 
each activity will be on the critical path. We shall call this proba- 
bility the Criticality Index , written, Computation of involves 

the concept of determining the probability that a particular p^* will be 
the maximum of a set of n random variables* j"p/*| . as before, these 
are not independent nor Identically distributed t hence this computation 
is not feasible analytically. It is not, however, difficult to deter- 
mine P c by the Monte Carlo method. This method also readily yields 
Tg, T l , and Fjy without resorting to the assumptions which were so 
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troublesome, but necessary, in the analytic approach, A detailed dis- 



cussion of this method, and suggestions for specific use of the criti- 
cality index will form a substantial portion of this study* 

The PERT Beta probability distribution was chosen as the distribu- 
tion for t& for an activity. Recall that in specifying this distribution 



By this assumption, we have specified the degree of uncertainty with 
which the estimator makes his three estimates. In other words* we are 
denying the difference between the case where the estimator may predict, 
with a high degree of confidence, that t will occur within a few time 
units of M s and the case where t has a higher likelihood of occurring 
near the extremes of the range of t*. The probability distribution for 
the former case should have a shape similar to that shown in Figure 5> 
while the latter case would be better represented by Figure 6„ It is 
not unreasonable to assume that, under certain conditions supervisors 
will have a basis for estimates with a narrow uncertainty band about the 
mode, even though the range of possible values of t* might be quite 
large* Under other conditions, the supervisor might be unable to place 
M with an appreciable degree of confidence In the interest of more 
accuracy in predicting the system behavior*, perhaps a choice of density 
functions should be available to the estimator to help describe his con- 
fidence in his estimates. This would not require a knowledge of proba- 
bility theory on the part of the supervisor. Hie could simply be asked 
to rate his uncertainty as to the relative position of M on a numerical 
scale, say with three descriptive choices* His choice would lead to a 
corresponding probability distribution. Appendices I, 12, HI and IV 



the variance of the distribution was arbitrarily chosen to be 
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describe means of providing this flexibility in choree of density func- 



tions with the desired degree of variance. 




Figure 5 

A Density Function with Small Variance 




Figure 6 

A Density Function with Large Variance 
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5, The Monte Carlo Technique 

We have seen that the completion time s t** for any activity in a 
proposed project is a random variable. As the project is performed* 
each activity is completed in some actual time, t. Only after comple- 
tion of the entire project can actual values of t, slack, and elapsed 
time, T, for events and activities be fixed, and the critical path de- 
termined. 

The Monte Carlo method calculates average values of these param- 
eters by repeatedly simulating the performance of the entire network, 
and statistically averaging the results. The simulation of the project 
is accomplished as follows: 

A random value of t for each activity in the network is drawn from 
the probability distribution for t* for the activity. These random 
values of t, called realizations of are drawn in such a way that in 
repeated drawings from a specific probability distribution, the number 
of realizations falling within a segment of fixed length on the t axis 
is proportional to the height of the probability density curve over that 
segment. Each activity has its own unique probability distri bu ticn* de- 
termined by the three estimates * hence, for one simulation of the project 
a realization of t* drawn for an activity is a completion time which 
could have resulted from a normal performance of the activity. One simu- 
lation of the project consists of drawing a value of t for each activity 
of the network. Calculations of the network are then performed, using 
these values of t instead of the t f} used in the analytical method, T and 
K for each activity are calculated, and then the critical path for that 
one simulation is determined. For each activity a tally is incremented 
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if the activity was actually on the critical path Another tally is in- 
cremented for each scheduled event if the event was accomplished prior 
to its scheduled completion date. 

The above simulation is repeated many times and the values of t, T, 
and K for each activity are added to cumulative totals after each Simula 
tion of the project* After an adequate number of replications, these 
cumulative totals and tallies are divided by the number ©f replications , 
yielding the average values 9 t e , Tg» and Kg for each activity and the 
probability that the activity will be on the critical path. For all 
scheduled events , the probability that the event was completed prior to 
the scheduled date is also computed* The observed variance of these 
statistics is also available if d@sired 9 as well as an approximate proba- 
bility distribution for project duration* 

The method of sampling the probability distributions of the activity 
times for realizations of t* is discussed in detail in Appendices 1, II, 
III, and IV. Using the above techniques, it is not difficult to provide 
a choice of probability distributions for t* in order to account for 
varying degrees of uncertainty in estimates, as was discussed earlier* 

The computer codes designed to implement the basic PERT system may 
be modified, without great difficulty, to perform the Monte Carlo calcu- 
lation* The author recently modified the North American Aviation PERT 
Mark III code in this manner* This program, designed for the IBM 7094, 
will accommodate a network of 1500 activities without requiring external 
memory units. Tie] 

The Monte Carlo technique is not subject t© the pitfalls ©f inde- 
pendence, assumptions of normalcy due t© the Central Limit Theorem, and 
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assumption of a particular longest path* etc. , which plague the analytic 
computation. The inaccuracies inherent in the Monte Carle technique are- 
due to non- randomness in the random number generator 9 and perhaps* fail- 
ure to perform a sufficient number of replications to permit the law of 
large numbers to obtain. Excellent random number generators are avail- 
able, providing results well within the limit ©f accuracy imposed by the 
time estimates. Statistical methods exist by which to determine the 
number of replications required in order to obtain the desired degree of 
confidence in the results. [l7^j 

The system is limited in accuracy only by the human error in estimat- 
ing activity times, and in approximating the correct probability distri- 
bution for t*. 

The Monte Carlo technique requires more computer running time than 
does the analytic method. An efficiently coded routine requires approxi- 
mately 750 microseconds to draw a random number and perform the calcula- 
tions necessary to transform it into a realization ©f t*. For details 
of the method employed, see Appendix 1. The approximate increase in 
PERT Program time required, if the Monte Carlo technique is used in lieu 
of the standard analytic method, may be estimated by the product 0 

(No. of activitxes)(No. of replic&t ions)(7S0 it ICT 6 sec). 

For example, a network consisting of 1000 activities would require 
12.5 minutes of additional computer running time if a Monte Carlo calcu- 
lation of 1000 replications were used instead of the standard technique. 
The advantages inherent in this technique would seem to render tnis in- 
crease in running time negligible * 
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6. The Extension to PERT/COST 10 

An extension of the basic PERT system , PERT/COST, was adopted 1 July 
1962 as the basic standard for management time and cost control systems 
by DOD and NASA. PERT/COST extends the original PERT system to provide 
cost estimation and cost monitoring features in addition to the procedure 
for prediction and control of time performance . Two additional optional 
features, a Time-Cost supplement , and a Resource Allocation supplement 
were also promulgated. 

PERT/COST adds to the basic system the following capabilities* 

a. Initial estimation of project cost, broken down by work pack* 
ages, sub-systems, etc. 

b. Provision of a consolidated estimate of the requirements for 
manpower and other resources as an aid to scheduling and pro- 
curement. 

c Provision of in progress cost reports, showing current cost 
status and revised cost estimates in comparison with budgeted 
expenditures and contract estimates. 

The optional Tirae«Cost supplement provides a method of estimating 
three feasible schedules for project completion, with their estimated 
costs and an evaluation of the technical risks associated with each 
schedule. These schedules are designed to provide completion in accord- 
ance with three criteria* (1) the most efficient use of time and re- 
sources i (2) completion by the DCD; (3) earliest feasible completion. 

The Resource Allocation supplement outlines a procedure for sched- 
uling a project to meet a DCD in an optimal manner with respect to costs, 
and utilization of manpower. 
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The concept of operation of PERT/COST is essentially as follows 
A* F r °. i cct breakdown . The project is broken down unto activities 
and work packages, and graphically represented by a network. For the 
purposes of cost reporting, estimating, and accounting, several activi- 
ties may be combined into a "work package 11 9 when the activities involved 
are small, and detailed cost reporting and estimating of these activities 
would involve unnecessary expense. In general, a work package does not 
comprise a closed network. That is, it may not be integrated and repre- 
sented by a single activity due to the presence of events interior to the 
work package which are connected by activities to other work packages. 

In this case, the work package must be broken down into its component 
activities and events for the network time calculations.. The cost esti- 
mating and expenditure reporting section of the system ©sy utilize the 
entire work package as the smallest organizational division. 

Certain work packages may not appear on the project network* In 
general, these are the functions of management not directed toward the 
accomplishment of specific activities and events. For example, the 
accounting, purchasing, or management travel costs for the project would 
be difficult to represent as activities in the network; however, the 
cost of these work packages may be estimated and reported,, and treated 
in the cost analysis portion of the PERT/COST system* 

B, Time c alc ul ations . Time estimates are made for each activity, 
and the network calculations performed. The project duration is com- 
pared with the DCDj, and if necessary* any of the following actions taken 
in order to reduce the overall project durations 

(a) Reorganize the network by increasing the parallelism of activi- 



ties. 
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(b) Commit additional resources to selected activities, or reallo- 
cate resources within the network. 

(c) Change or delete activities. 

C. Manpower loading . [9] A breakdown by skills of manpower required 
for an activity/work package may be included with the time and cost esti- 
mates e Computer programs now in existence assemble this information and 
prepare a report and display showing the project manpower requirements, 

by skills, plotted against time. 7, llj This information provides project 
management with a forecast of manpower needs, and a basis for reschedul- 
ing activities in order to best utilize available manpower. Scheduling 
a large project without considering the overall use of manpower skills 
leads to very uneconomical manpower use. One time period may require a 
manpower level for a certain skill far in excess of that available, while 
the following time period makes little use of that skill in performing 
the necessary rescheduling, the time constraints of the network and 
criticalness of activities must be considered. Those activities least 
likely to be critical are the obvious candidates for rescheduling for 
the purpose of leveling manpower requirements. Computer programs now in 
existence, supplementing standard PERT routines, perform this reschedul- 
ing under the constraints of maximum availability of manpower and with 
consideration of slack. 7, llj 

D. Scheduling the network . A scheduled duration, t s * is promul- 
gated for each activity. These times may be less than, equal to, or 
greater than t e , depending on the situation and management policy. Using 
the t s for each activity, the network is computed and values of 3g and 
S^, corresponding to Tg and Tg, determined. Sg represents the earliest 
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date on which an activity/event can be scheduled for completion, 
represents the latest date on which an activity/event may be scheduled 
for completion without causing a schedule slippage for the project com- 
pletion. and S L are computed in exactly the same manner as were 

and when t g is used as activity time instead of t 0 . As a final step 
in scheduling, an elapsed time, T s , is imposed at each event, T s is 
chosen to lie between Sg and S^, and may be translated late a firm sched 
uled completion date for activities immediately preceding the event, and 
a beginning date for activities immediately following the event, 

E. Budgeting the project . After the schedule is finalized, cost 
estimates are prepared for each work package , and a budget prepared for 
allocating funds, 

F, Koflltoring progress of the project . As the project progresses s 
reports are submitted from activities/work packages 8 indicating progress 
new time estimates, funds committed to date, and new estimates for cost. 
This information forms inputs to the PERT/COST computer program, which 
performs the network calculations, and prepares reports showing new 
values for Tg, T^, Kg, P D , and the new critical path s as well as current 
expenditures in comparison with estimates and budget. The new projected 
cost curve is also generated and compared with estimates and budget. 
Armed with this information, management may reschedule, reallocates, re- 
budget, or take whatever action as may be appropriate to the situation. 

The Time -Cos t Option Supplement . The purpose of this supplement is 
to standardize procedures by which a project manager may prepare three 
alternative project schedule proposals, and evaluate the technical risks 
associated with each. By technical risk, we mean the gamble of 
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performance, time, and cost incurred by departing from the best develop- 
ment technique in order to meet the schedule * 

First a plan is formulated designed to meet the project DCD, as 
imposed by the controlling agency. 

Next, a Most Efficient plan is formulated, in which the best devel- 
opment techniques are used in order to reduce technical risks, and most 
efficient use is made of manpower and resources, resulting in greater 
duration for the project. Total cost is generally lower for this plan. 

Finally, a Shortest Time plan is proposed. The purpose of such a 
plan is to explore the feasibility of scheduling the project for a dura- 
tion shorter than that proposed by the contracting agency in order to 
benefit from any strategic benefits which might accrue from such early 
completion. In order to arrive at the shortest time plan, the schedule 
is compressed by allocation of additional resources, paralleling activi- 
ties which should normally go in sequence, eliminating activities, and 
changing technical approaches to the problems of design and construction. 
All such actions serve to increase the technical risks involved, and most 
of them increase cost. 

The three options are presented to the contracting agency with esti- 
mates of cost and time, and with a comparison and evaluation of the tech- 
nical risks associated with each plan. 

The PERT/COST system is the primary tool for obtaining the time and 
cost estimates needed to prepare the three schedule proposals. 

The Resource Allocation Supplement , This supplement provides a pro- 
cedure by which project managers may schedule the project in the optimal 
manner with regard to the time-cost trade-off. The concept of operation 
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is as follows. 



For any activity there may exist several feasible schedules, mvolv 
ing different time-cost relationships „ (See Figure 7,) Efficient use 
of manpower and existing resources and machinery may result in schedule 
A. Further time extension to B might result in higher costs due to the 
effect of fixed costs. By use of overtime, hiring additional manpower, 
or purchase of additional machinery or space, time-cost combination C or 
D might be achieved. 

Cost h 
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Figure 7 

Time-Cost Plans for an Activity 

Schedule A would be the most economical for this activity, hence 
would be the time estimated for the activity for the first network cal- 
culation, in a similar manner, most economical estimates are obtained, 
with corresponding cost estimates, from the other activities „ The net- 
work calculation is then performed, and the project duration,, D, compared 
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with Dg . If D>Dg, select ed activities on the critical path must be 
shortened by resorting to a higher cost-shorter time work schedule 

For each activity on the critical path, determine the increase in 
cost, divided by the decrease in time that would occur in moving to eac b 
shorter time-cost point. This is the slope of the time cost curve for 
the corresponding time reduction. The activity with the lowest value of 
slope is then chosen as the activity to reschedule. 

If this time reduction in the critical path exceeds the difference 
in slack between the critical path and the smallest slack not on the 
critical path, it becomes necessary to recompute the network and deter- 
mine a new critical path. The above process is repeated as often as 
necessary until D ^ Dg. 

Slack paths are then reexamined to see if any activities may be ex- 
tended to lower cost points without going critical. If fixed costs were 
not a part of the activity estimates, notice may be taken of the fact th 
the most economical point to operate is usually to the left of the mini- 
mum cost point with respect only to direct costs. After determining the 
most economical plan with respect to direct costs, a test is made to de- 
termine if further time reduction below the DCD would result in lower 
overall costs. This is done bv adding on the fixed costs for the project 
and attempting further reduction in time. If further time reduction re- 
sults in lower overall cost, then further time reduction is performed 
until total costs reach a minimum. 

The schedule is then adjusted as necessary to reach an optimum lev- 
eling of manpower and other resources within the existing constraints 
and considering total cost. 
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7. A Method for Resource Allocation, using the Criticality Index 

In Section 4, we pointed out the errors in the concept of a deter- 
ministic prediction of the critical path. We also showed the optimistic 
bias inherent in the assumption that the longest path to an event will 
be the path whose mean length is greatest. The Monte Carlo technique 
does not depend on this assumption, hence, with this technique we can 
calculate Pq, Tg, and with accuracy limited only by the supervisor's 
estimates and our approximation of the true probability distributions of 
t* for the activities. We shall now develop a system for optimum resource 
allocation utilizing the criticality index obtained from the Monte Carlo 
calculation. 

In order to use this method, we must require that a work package be 
a glo se .d network . A closed network is a network with a single beginning 
event and a single end event. Activities exterior to the closed network 
may not be incident to, or emanate from, events within the closed network 
other than the beginning or end event. (See Figure 8.) In this example, 
if any activities not included in the work package should be incident to, 
or emanate from events B, C, or D, then the network would not be closed. 
This same restriction is implicit in the resource allocation supplement 
of the standard PERT/COST system, although not stated in the PERT/COST 
document. Notice that the closed network may be integrated and repre- 
sented by a single activity connecting events A and E. Some PERT com- 
puter programs now in existence are capable of integrating such sub- 
networks for simplification of the network for the benefit of higher 
levels of management. JJLSj 
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Figure 8 

A Closed Network 

We next require that the supervisor in charge of each work package 
submit as many time-cost plans as are feasible for his work package com- 
pletion, These plans are represented by points on the time-cost curve in 
the same manner as previously discussed in connection with the Resource 
allocation supplement. No points to the right of the minimum cost point 
are submitted. The time-cost information is represented as a set of n 
pairs, (time, cost), in descending order of time and corresponding ascend- 
ing order of cost. For example: 

14 weeks , $40,000 
12 weeks, 60,000 
10 weeks , 90,000 

The time corresponding to the minimum cost point is taken as the 
most likely estimate, M, of activity completion time. The additional 
time estimates, $ and P, are also submitted for the minimum cost plan. 
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These three estimates and the corresponding cost are the basis of the 
preliminary calculation of the network. From this calculation, the prob- 
ability, P^, of meeting the project DCD, is determined. Recall that this 
probability is considerably more valid as a result of the Monte Carlo 
technique . 

If Pj) is unacceptably small, the network must be compressed by re- 
designing the network layout and/or allocating more resources to selected 
work packages in order to achieve shorter completion time. If, after 
whatever network redesign is considered feasible, the network must be 
further compressed, the resource allocation phase of the program is per- 
formed. This phase consists of the following logical steps: 

1. For each time-cost plan of each activity having positive P q v 
the following figure of merit is computed. 

etb = p c 2^3 (*) 

where /\t is the reduction in time achieved by using the T-C plan being 
considered, over that used in the previous network calculation. AC is 
the corresponding increase in cost. ETB may be considered as the expected 
amount of time which will be "bought" per unit cost by selecting the T-C 
plan being considered. 

2 . Form 3-tuples consisting of each ETB, its corresponding At, 
and its activity identification. Order this list in descending order of 
ETB. 

3. Select from the top of the list, a sufficient number of reduced 
time plans so that the sura of these t's does not exceed a predetermined 
percentage, p, of D. The best value for p may be determined experimen- 
tally. The selection is made from the top of the list, but only one 
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plan, that with the largest At, is retained for any one activity, 

4. For the activities chosen for time reduction, the time corre- 
sponding to the new T-C plan is taken as the new value of M, The esti- 
mates for P and 0 are revised downward by multiplying by the ratio, 

. (This procedure may be subject to question; another method 

M (old) 

would be to require P and 0 estimates for each T-C plan,) 

5, With the new time and cost estimates for the accelerated activi- 
ties, the network is again calculated, and P|) determined. All activities 
will have a new P^, The mean duration, and P D are again evaluated, if 
still unacceptable, the process may be continued until either an accept- 
able Pq is attained, or the budgetary limit is reached. The schedule is 
finalized on the basis of the time plans used in the last network calcu- 
lation, after consideration of manpower and other resource leveling. 

By the above method, resources are allocated where, probabilistic- 
ally, they can be expected to contribute most effectively to shortening 
project duration. The requirement for several feasible time-cost plans 
may serve the additional purpose of forcing line supervisory personnel 
to consider the time-cost relationship more carefully, allowing operation 
at the most efficient point on the time-cost curve whenever feasible. 

An optimum value for p can be determined with a few experimental 
runs of the program, p should be chosen as large as possible in order 
to decrease the number of network calculations required, but small enough 
so that successive calculations produce a smooth compression of the net- 
work. Operation of the program will tend to decrease the differences 
between the ETB*s for the activities, which is equivalent to operating 
the activities at the same level of expected marginal utility. This is 
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a well known principle of optimization* 

The inputs to the system, three time estimates and a T-C function 
specified by several 2- tuples, are not complicated, and are compatible 
with the concept of progress and cost monitoring envisioned in the PERT/ 
COST system. 

A flow diagram of the logic used in this resource allocation method 
is shown in Appendix V. 
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8. A Military Application of the PERT Network Approach 

One of the most outstanding features of the PERT approach is the 
graphical portrayal of the project by means of the network. By this de- 
vice, relationships between the various tasks may be clearly visualized. 
The insight thus gained by the managers of the project, into these com- 
plex relationships, enables them to function with greatly increased 
effectiveness in their managerial capacity. 

A military operation is in many ways similar to a development proj- 
ect. The operation may be divided into a number of tasks, with many 
interrelationships existing between individual tasks. Hence, the net- 
work may serve the same useful function in this application as in indus- 
trial development. 

The phases of military planning known as "The Development @f the 
Plan," and "Supervision of the Planned Action," are most readily bene- 
fited by the network approach. In the formalized military planning 
sequence, the Development of the Plan phase has been preceded by the 
"Estimate of the Situation," in which the assigned mission has been 
studied, possible outcomes based on enemy capabilities and alternative 
own courses of action have been analyzed in a game theoretic matrix, 
and a decision has been reached regarding a general course of action to 
be pursued. 

The problem now involves developing a complete, detailed plan for 
the operation, including organizing the available forces and assigning to 
each task unit the appropriate tasks which make up the general plan. At 
this stage, detailed planning is done in which training, acquisition of 
intelligence, movement, communications, logistics, and battle action are 
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all considered. 



The plan begins with an initial concept, which may be represented 
by a relatively simple network, whose tasks are stated in broad terms of 
accomplishment. As the planning becomes more detailed, these initial 
tasks may be subdivided into networks of lower-order tasks. This ph s 
may be done at lower echelons of command, upon receipt by them of their 
superior 1 s directive . 

The network will aid in establishing a firm schedule for the opera 
tion, in much the same manner as is done in industrial projects. It: aay 
also be a valuable aid for conducting briefings for subordinates, since 
it will enable them to better appreciate the relationship of their own 
tasks to the whole operation,, 

During the operation, the network may serve as a display, on which 
current progress is indicated as reported. The effects of delays and 
failures may be more readily evaluated and corrected* 

The initial concept of a plan to execute an amphibious landing in 
eneray-held territory might be represented by the network show, in Figure 9 
A network may also be used to lay out the basic strategy of a pro- 
tracted military campaign. In this usage, tasks would be represented by 
subsidiary operations, such as: "Seize Island B, ,T or "Cut supply line 

XY." The procedure is similar to that employed in the operation just dis- 
cussed, except that it takes place over an extended period of time* In 
many cases, proposed actions are tentative, or the selection of alterna- 
tive actions may depend on future developments, such as the success or 
failure of a previous endeavor, or on the subsequent enemy reaction. We 
may indicate this need for decision at a particular point in the network 
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by showing a "decision box" in lieu of an event at that point. The deci- 
sion box has more than one task emanating from it; the decision, to be 
made at some later date, determining which task(s) to be performed. The 
network then proceeds to show the planned action to be taken for each 
alternative. Figure 10 is an example of a portion of such a strategy 
network. [l3,14*[ 
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APPENDIX I 



Monte Carlo Techniques with the Beta Distribution 
The model for activity completion time, t*, chosen by the originators 
of PERT, is the Beta probability distribution, characterized by its den- 
sity function: [ 2j 

£ t*(t) = K(t-0)^ (P-t/ , 0 - t — P (5) 

= 0, elsewhere. 

We may reduce this to the standard form of the Beta distribution by the 

following transformation: 

x * = t*- 0 
P - 0 



The probability density function of x* is: 

f x *(x) = Kj_ x x (1-x) , 0 — t — 1 



( 7 ) 



The variable, x*, is equivalent to t*, where 0=0 and P = 1. 

The two parameters, 0 and P, specify the upper and lower limits of 

the variation of t*. It is now desired to specify the remaining two 

parameters , and / , so that the mode of the density function occurs at 

M, and so that the distribution has the desired variance. 

Let ra « denote the mode of x*. By finding the root of f* (x), 

P-0 x 



we obtain: 



m = 



c<+ f 

The variance of x* is given by: 
2 



( 8 ) 



, 2 (<*•+ 1) (f+ 1) 

0 ^ X ‘'^"(oC+iT'-+-2) 2 (oC+^'+3) 



(9) 



Now suppose we desire that the standard deviation, cf (t*), he expressed 

/p . ^ \ j \ 

by/— I, where d is an arbitrary constant. Then Q (x*) = — . By 
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requiring that 0 ( x* ) - “f ’ we ohta ^ n t ^ ie following relationship in c >4 



and in 
3 

OL 4 



y 2 in 3 .dV b 7m) oc* + - J 



( 10 ) 



Now ©i can be determined as a function of m and d by the above equation 
With o( determined, I ( is specified by solving equation ( 8 ): 



m 



The mean of x* is given by; 

oi + 1 



E[xj . 



oC + f +£. 



( 11 ) 



( 12 ) 



Since must be determined from the cubic equation, (10), for each value 
of m, the computation of E(Vj is cumbersome. In the PERT Beta model, d 
was chosen as six. For this value of d, E[xJ is approximately linear in 
m, and may be approximated by the following relation: 



E[x] • $”+_! . 

o 

The transformation to t* yields: 

0 ± 4M + P 
6 




(13) 



(14) 



Random Sampling from the Beta Distribution . There are several good 
random number generating routines which yield random numbers from the uni- 
form (0,1) distribution. We shall use the technique of the Probability 
Integral Transformation to transform the random number, u, drawn from the 
uniform (0,1) distribution, to a corresponding sample from whatever other 
probability distribution we desire. [ 3 ] 

This transformation depends on the following theorem; 
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Fur any ran dom v ariable , v* , hav ing JLhe prob abili ty distri bution 
fun ction , F v *, and the density f unction , defi ne the r andom variable , 

u* - F that is: 



u 



F v *< v) s 



f .. (C)dt 
V* 



/ ^.oO 



( 15 ) 



then , u* ijs a uniformly distributed random variable on the interval (0,1) 
Hence, we may use the inverse function, v = F ^(u), to transform the 
sample, u, from the uniform (0,1) distribution to the corresponding sample , 
v, from the desired distribution. 

The distribution function: 



F x *< x > 




oL tf 
y (l-y) dy 



( 16 ) 



for the Beta variable, x*, must be calculated by numerical integration 
In order to utilize the procedure given above for sampling from this dis- 
tribution, we may store the tabled distribution functions in the memory 
of the computer and use a table look-up routine to enter the table with 
the random sample, u, and determine the corresponding value for x 0 A 
complete distribution function must be stored for each of the increments 
of m, and for a limited number of values of d. 

The input card for an activity specifies the choice of d which best 
represents the estimator's uncertainty. Also specified are the three 
estimates, $, M, and P, With d and M specified, the proper table may be 
selected, u is produced by the random number generator. The correspond- 
ing value of x is obtained from the table, then transformed into a reali- 
zation of t* by the following transformation: 
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t = (P-0) X -f 0 



(17) 

Tables 1 through V are tables of distribution functions for tr t 

standard Beta distribution, with d taking on values of 4 through 8, re 

spectively. A complete function is given for each of 11 values of m 

As an example of the use of the procedure given above, suppose 0, 

and P were 12, 15, and 20 weeks, respectively, and the selection, d * 6 V 

15-12 

is made. Then ra =• — r— = *375. Suppose we draw the random number, 

o 

.795, from the uniform random number generator. Entering Table III for 
m = .3, with the argument, .795, we obtain x = ,50. From the table for 
os .4, linear interpolation yields x = .572. Interpolating between 
these two values of x, we obtain x = .554 for m = .375. The correspond- 
ing realization of t* is then 8 x -f 12 r 16.4. 

An efficiently coded routine designed to perform the above table 
searches and calculations, requires approximately 750 microseconds. 

Figures 11, 12, and 13 are families of curves of the probability 
density functions for x* corresponding to d r 4, 6, and 8, respectively. 
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DISTRIBUTION FUNCTION OF THE STANDARD BETA FUNCTION WITH D = 





A= .000 


A= .079 


A= .174 


A= 


.280 


A= .393 


A= . 500 




G= .720 


G= .714 


G= .694 


G= 


.654 


G= .589 


G= . 500 


X 


M= .000 


M= .100 


M= .200 


M= 


. 300 


M= .400 


M= .500 


.05 


.071 


.059 


.047 




.035 


.025 


. 018 


• 10 


.068 


.062 


.054 




.045 


.037 


. 029 


. 1 5 


.233 


.207 


. 179 




. 1 48 


. 1 19 


. 093 


.20 


.309 


.280 


.247 




.211 


.175 


. 142 


.25 


.381 


.352 


.316 




.2 76 


.234 


. 195 


• 3 0 


.451 


.421 


.384 




. 342 


.296 


.252 


• 35 


.516 


.487 


.45 1 




. 407 


.360 


.312 


• 40 


.579 


.551 


.515 




.473 


.424 


. 373 


• 45 


. 6 37 


.611 


.578 




.536 


.488 


. 436 


• 50 


.692 


.668 


.638 




. 598 


.552 


. 500 


• 55 


.743 


.722 


.694 




. 658 


.614 


. 564 


• 60 


.790 


‘ .772 


.748 




.715 


.674 


. 627 


• 65 


.833 


.818 


.797 




.769 


.732 


. 688 


• 70 


.872 


.860 


.843 




.819 


.787 


. 748 


• 75 


.907 


.897 


.884 




. 865 


.839 


. 805 


• 00 


.936 


.930 


.920 




. 906 


.885 


. 858 


• 85 


.961 


.957 


.951 




.941 


.927 


. 907 


• 90 


.981 


.979 


.975 




.970 


.961 


. 949 


• 95 


.994 


.994 


.992 




.990 


.987 


. 982 


1*00 


1.000 


1.000 


1.000 




1.000 


1.000 


1. 000 





A= .590 


A= .654 


A= .694 


A= 


.714 


A= .720 




G= .393 


G= .280 


G= .174 


G= 


. 079 


G= .000 


X 


M= .600 


M= .700 


M= .800 


y= 


. 900 


M= 1 . 000 


.05 


.013 


.010 


.008 




. 006 


.006 


.10 


.023 


.018 


.015 




.013 


.012 


. 1 5 


.073 


.059 


.049 




. 043 


.038 


.20 


.115 


.094 


.080 




. 070 


.063 


.25 


. 161 


.135 


.116 




. 102 


.092 


.30 


.213 


. 181 


.157 




. 1 40 


. 1 26 


.35 


.268 


.231 


.203 




.181 


.164 


.40 


.326 


.285 


.252 




.227 


.207 


.45 


.386 


.342 


.306 




.277 


.253 


.50 


.448 


.402 


.362 




. 331 


.304 


.55 


.512 


.464 


.422 




. 388 


.358 


.60 


.576 


.527 


.484 




.4.48 


.415 


. 65 


.640 


.593 


. 549 




.512 


.477 


.70 


.704 


.658 


.616 




. 578 


.541 


.75 


.766 


.724 


.684 




. 647 


.610 


.80 


.825 


.789 


.753 




.718 


.681 


.85 


.881 . 


.852 


.821 




. 791 


.756 


.90 


.932 


.911 


.889 




.865 


.834 


.95 


.975 


.965 


.953 




.938 


.916 


1.00 


1.000 


1.000 


1.000 




1.000 


1.000 



^Table I 

. f/ U-y) G ay 

F x .(x) - ~~~7TT r 

x d-y) G ay 



M is the position of the mode 
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01 ST R 1 8UT ION FUNCTION OF THE STANDARD BETA FUNCTION WITH 0 = 5 





A= .000 


A= .217 


A= .509 


A= .870 


A= 1 .262 


A=1. 625 




G=1 .826 


G=1 .951 


G=2.036 


G=2. 030 


G=1 .893 


G*l. 625 


X 


M= .000 


M= . 100 


M= .200 


M= .300 


M= .400 


M* . 500 


.05 


.114 


.078 


.044 


.019 


.007 


. 002 


.10 


.104 


.088 


.066 


. 041 


.022 


.010 


.15 


.353 


.288 


.210 


. 132 


.072 


. 036 


.20 


.455 


. 389 


‘ .304 


.210 


.129 


' . 072 


.25 


.546 


.484 


.398 


.296 


.197 


. 120 


.30 


.626 


.570 


.489 


. 385 


.275 


. 181 


.35 


.697 


.648 


.575 


.475 


.360 


. 251 


.40 


.758 


.718 


.653 


. 561 


.447 


. 329 


.45 


.811 


.778 


.724 


.643 


.534 


.413 


.50 


.856 


.830 


.7e7 


.717 


.619 


. 500 


.55 


.893 


. 874 


.840 


.784 


.698 


. 587 


.60 


.923 


- .910 


.885 


.841 


.770 


. 671 


.65 


.947 


.939 


.921 


. 889 


.833 


. 749 


.70 


.966 


.961 


.949 


. 927 


.887 


.819 


.75 


.980 


.977 


.970 


. 956 


.929 


. 880 


.80 


.989 


.988 


.985 


. 977 


.961 


. 928 


.85 


.995 


.995 


.993 


. 990 


.982 


. 964 


.90 


.998 


.998 


.998 


. 997 


.994 


. 987 


.95 


1.000 


1.000 


1.000 


1. 000 


.999 


. 998 


1 .00 


1.000 


1.000 


1.000 


1. 000 


1.000 


1. 000 




A= 1 . 894 


% 

A=2 • 030 


A=2. 036 


A= 1. 9 5 1 


A= 1 • 826 






G=1 .263 


G= .870 


G= .509 


G= .217 


G= .000 




X 


M= .600 


M= .700 


M= .800 


K= .900 


M=1 .000 




.05 


.001 


.000 


.000 


. 000 


.000 




.10 


.005 


.002 


.002 


.001 


.001 




.15 


.018 


.010 


.007 


. 005 


.005 




.20 


.039 


.023 


.015 


.012 


.011 




.25 


.071 


.044 


.030 


. 023 


.020 




.30 


.113 


.073 


.051 


.039 


.033 




.35 


. 166 


.111 


.079 


. 06 1 


.051 




.40 


.230 


.159 


.115 


. 090 


.075 




.45 


.302 


.216 


.160 


. 126 


.105 




.50 


.381 


.283 


.213 


. 170 


.141 




.55 


• 466 


.357 


.276 


. 222 


.185 




.60 


.553 


.439 


.347 


. 282 


.236 




.65 


.640 


.525 


.425 


. 352 


.296 




.70 


.725 


.615 


.511 


.430 


.365 




.75 


.803 


.704 


.602 


.516 


.444 




.80 


.871 


.790 


.696 


.611 


.532 




.85 


.928 


.868 


.790 


.712 


.632 




.90 


.969 


.934 


.879 


.817 


.742 




.95 


.993 


.981 


.956 


.922 


.865 




1.00 


1.000 


1.000 


1.000 


1.000 


1.000 








Table 

V 


II 










II 

/-N 

X 

* 

X 


1/ a- 


,y) G dy 






















f 1 y A U-y) G dy 









'o 



M is the position of the mode 
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, **^9/%a* 






DISTRIBUTION FUNCTION OF THE STANDARD BETA FUNCTION WITH 0=6 





A= .000 


A= .361 


A= .891 


A= 1. 574 


A=2.320 


A=3. 000 




G=2 . 072 


G=3. 252 


G=3. 565 


G=3. 673 


G=3. 480 


G=3. 000 


X 


M= .000 


M= . 100 


Ha .200 


M= .300 


M= .400 


M= . 500 


.05 


.153 


.090 


.035 


. 008 


.001 


. COO 


.10 


.131 


. 109 


.071 


.032 


.010 


. 002 


.15 


.449 


.350 


.222 


. 107 


.040 


.012 


.20 


.564 


.472 


.338 


. 195 


.088 


. 033 


.25 


.661 


.581 


.454 


. 298 


.159 


.071 


.30 


.740 


.675 


.563 


.408 


.248 


. 126 


.35 


.805 


.755 


.662 


.5 19 


.350 


. 200 


.40 


.857 


.820 


.747 


.624 


.458 


. 290 


.45 


.89e 


.072 


.818 


.718 


.567 


.392 


.50 


.929 


.912 


.874 


.798 


.670 


. 500 


.55 


.953 


.943 . 


.917 


.863 


.761 


. 608 


.60 


.970 


- .964 


.949 


.913 


.838 


.710 


.65 


.982 


.979 


.971 


. 949 


.898 


. 800 


.70 


.990 


.989 


.985 


.973 


.942 


. 874 


.75 


.995 


.995 


.993 


. 987 


.971 


. 929 


.80 


.998 


.998 


.997 


. 995 


.988 


. 967 


.85 


.999 


.999 


.999 


. 999 


.996 


.988 


.90 


1.000 


1.000 


1.000 


1.000 


.999 


. 997 


.95 


1.000 


1.000 


1.000 


1.000 


1.000 


1. 000 


1.00 


1.000 


1.000 


1.000 


1. 000 


1.000 


1. 000 





A=3.48 1 


A=3. 672 


A=3.565 


A=3* 2 52 


A=2. 872 




G=2 .321 


G=1 . 574 


G= .891 


G= .361 


G= .000 


X 


M= .600 


M= .700 


M= .800 


.900 


M= 1 . 000 


.05 


.000 


.000 


.000 


.000 


.000 


.10 


.001 


.000 


.000 


.000 


.000 


.15 


.004 


.001 


.001 


. 001 


.001 


.20 


.012 


.005 


.003 


. 002 


• .002 


.25 


.029 


.013 


.007 


.005 


.005 


.30 


.058 


.027 


.015 


.011 


.009 


.35 


.102 


.051 


.029 


.021 


.017 


.40 


.162 


.087 


.051 


.036 


.029 


.45 


.239 


. 137 


.083 


.057 


.045 


.50 


.330 


.202 


.126 


.088 


.068 


.55 


.433 


.282 


.182 


. 128 


.099 


.60 


.542 


.376 


.253 


. 180 


.138 


.65 


.650 


.481 


.338 


.245 


.189 


.70 


.752 


.592 


.437 


.325 


.251 


.75 


.841 


.702 


.546 


.419 


.328 


.80 


.912 


.805 


.662 


. 528 


.421 


.85 


.960 


.893 ' 


.778 


.650 


.533 


.90 


.988 


.956 


.884 


.781 


.665 


.95 


.999 


.992 


.965 


.910 


.820 


1.00 


1 .000 


1.000 


1.000 


1.000 


1.000 



F x+ (x) 



Table III 




a y 
a y 



M is the position of the mode. 
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DISTRIBUTION FUNCTION OF THE STANDARD BETA FUNCTION WITH D * 

A = .000 
G=3 . 899 
X M= .000 



1 



.05 
. 10 
.15 
.20 
.25 
.30 
.35 
.40 
.45 
.50 
.55 
.60 
.65 
.70 
.75 
.80 
.85 
.90 
.95 
.00 



.189 
.154 
.530 
.651 
.745 
.818 
.874 
.915 
.944 
.965 
.979 
.988 
.994 
.997 
.999 
1.000 
1 .000 
1.000 
1.000 
1.000 



A= .517 


A* 1.329 


A=2. 399 


A=3. 569 


A=4. 625 


G=4. 655 


G=5 . 3 1 6 


G=5. 598 


G=5.353 


G=4. 625 


M= .100 


M= .200 


.300 


M= .400 


M= . 500 


.097 


.026 


. 003 


.000 


. 000 


. 128 


.071 


.022 


.004 


. 000 


.401 


.225 


. 083 


.020 


. 004 


. 540 


.361 


. 175 


.058 


.014 


.658 


.497 


.292 


.124 


. 039 


.754 


.622 


.423 


.219 


. 085 


.830 


.729 


. 554 


.336 


. 157 


.886 


.815 


. 675 


.465 


.254 


.927 


.881 


. 777 


.595 


. 371 


* 956 


.928 


.857 


.713 


. 500 


.975 


.959 


.915 


.812 


. 629 


- .987 


.979 


.954 


.888 


. 746 


.993 


.990 


. 978 


.940 


. 843 


.997 


.996 


. 991 


.972 


.915 


.999 


.999 


. 997 


.989 


. 961 


1 . 000 


1 .000 


. 999 


.997 


. 986 


1.000 


1.000 


1.000 


.999 


. 996 


1.000 


1.000 


1. 000 


1.000 


1. 000 


1 . 000 


1.000 


1. 000 


1.000 


1 . 000 


1.000 


1.000 


1. 000 


1.000 


1. 000 



.05 

.10 

.15 

.20 

.25 

.30 

.35 

.40 

.45 

.50 

.55 

.60 

.65 

.70 

.75 

.80 

.85 

.90 

.95 

1.00 



A=5 . 354 


A=5. 597 


A=5. 3 1 6 


A=4. 656 


G=3. 569 


G=2 . 399 


G= 1.329 


G= .517 


M= ..600 


M= .700 


M= . 8C0 


M= .900 


.000 


.000 


.000 


. 000 


.000 


.000 


.000 


. 000 


.001 


.000 


.000 


. 000 


.003 


.001 


.000 


. 000 


.011 


.003 


.001 


.001 


.028 


.009 


.004 


. 003 


.060 


.022 


.010 


. 007 


.112 


.046 


.021 


.013 


.188 


.085 


.041 


.025 


.287 


. 143 


.072 


. 044 


.405 


.223 


. 119 


. 073 


.535 


.325 


.185 


.114 


.664 


.446 


.271 


. 170 


.78 1 


.577 


.378 


.246 


.876 


.708 


.503 


. 342 


.942 


.825 


.639 


.460 


.980 


.917 


.775 


.599 


.996 


.974 


.894 


. 752 


1.000 


.997 


.974 


. 903 


1.000 


1.000 


1.000 ' 


1. 000 



A=3. 899 
G= .000 
M= 1.000 

.000 

.000 

.000 

..000 

.001 

.003 

.006 

.011 

.020 

.034 

.053 

.082 

.121 

.174 

.244 

.335 

.451 

.597 

.778 

1.000 



Table IV 



F x „(x) 



J/ (1-7) 5 d 7 
(i-y) G dy 



M is the position of the mode 
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DISTRIBUTION FUNCTION OF THE STANDARO BETA FUNCTION WITH D = 8 





A* .000 


A= .686 


A= 1.825 


A=3. 348 


A=5.010 


A=6. 050 




G=4 .916 


G=6. 175 


G=7. 300 


G=7. 8 1 2 


G=7. 5 1 5 


G=6. 050 


X 


M= .000 


M= .100 


M= .200 


H s . 300 


M* .400 


M= .500 


.05 


.223 


. 101 


.0 18 


.001 


.000 


. 000 


.10 


. 1 12 


.145 


.069 


.013 


.001 


. 000 


.15 


.598 


.445 


.223 


. 063 


.009 


. 001 


.20 


.719 


.598 


.378 


. 1 54 


.036 


. 007 


.25 


.808 


.721 


.534 


.283 


.095 


. 024 


.30 


.872 


.815 


.672 


.434 


.192 


. 062 


.35 


.918 


.883 


.784 


. 585 


.321 


. 129 


.40 


.949 


.929 


.867 


.718 


.470 


. 228 


.45 


.969 


.960 


.924 


.825 


.620 


. 356 


.50 


.903 


.979 


.960 


. 901 


.751 


. 500 


.55 


.991 


.989 


.981 


. 949 


.854 


. 644 


.60 


.995 


‘ .995 


.992 


.977 


.924 


. 772 


.65 


.998 


.998 


.997 


. 991 


.966 


. 871 


.70 


.999 


.999 


.999 


. 997 


.988 


. 938 


.75 


1.000 


1.000 


1.000 


.999 


.997 


. 976 


.00 


1.000 


1.000 


1.000 


1.000 


.999 


. 993 


.85 


1.000 


1.000 


1.000 


1.000 


'1.000 


. 999 


.90 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


.95 


1.000 


1.000 


1.000 


1. 000 


1.000 


1. 000 


1.00 


1 .000 


1.000 


1.000 


1.000 


1.000 


1. 000 





A=7 .511 


A»7 . 8 1 1 


A = 7. 0 30 


A = 6. 175 


A= 4. 9 1 6 




G=5. 007 


G=3 . 348 


G= 1 . 758 


G= .686 


G= .000 


X 


M= .600 


M= .700 


M* .800 


M= .900 


M* 1 .000 


.05 


.000 


.000 


.000 


.000 


.000 


.10 


.000 


.000 


.000 


. 000 


.000 


.15 


.000 


.000 


.000 


.000 


.000 


.20 


.001 


.000 


.000 


.000 


.000 


.25 


.003 


.001 


.000 


.000 


.000 


.30 


.012 


.003 


.001 


.001 


.001 


.35 


.034 


.009 


.004 


.002 


.002 


.40 


.076 


.023 


.009 


. 005 


.004 


.45 


.146 


.051 


.021 


.011 


.009 


.50 


.249 


.099 


.043 


.021 


.017 


.55 


.380 


. 175 


.081 


.040 


.029 


.60 


.530 


.282 


.139 


.071 


.049 


.65 


.679 


.415 


.222 


. 1 17 


.078 


.70 


.808 


.566 


.334 


. 185 


.121 


.75 


.905 


.717 


.470 


.279 


.182 


.80 


.964 


.846 


.624 


.402 


.267 


.85 


.991 


.937 


.777 


. 555 


.382 


.90 


.999 


.985 


.904 


. 728 


.536 


.95 


1.000 


.999 


.981 


. 899 


.738 


1.00 


.... 1.000 


1.000 


1.000 


1.000 


1.000 


< • 




’ >1 f. i • - > •» 


"" 1 




• 



Table V 



T? (v') - 


(i-y)° dy 


" r 


V d-y) G dy 



'a 



M is the position of the mode 
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Beta Density Functions for d * A 
Modes at 0 9 0 2 s and <>5 
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Figure 12 

Beta Density Functions for d = 6 
Modes at 0, . 1, .2, «3, .4, and .5 
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Figure 13 

Beta Density Functions 
Modes at 0 9 . l 9 ,2, .3, 



for d = 8 
„4, and «5 
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APPENDIX II 



A Simpler Method for Operating With the Beta Distribution 
Since the Monte Carlo method does not require that the variance of 
t* be calculated and applied in order to evaluate P^, we may consider dis- 
pensing with the requirement that the variance of t* be constant over the 
range of M. This requirement is not necessarily a natural one. The vari- 
ance of the triangular distribution, for example, defined naturally by 0, 
M, and P, is given by: 

cr 2 (t*) = (l8) 



We see chat the variance for this distribution varies from a maximum of 

with m at the extremes of its range, to — (AiA) with m at the 
midpoint. The variance of the SinCos function, (Appendix IV), another 
function naturally specified by the three estimates, behaves in the same 
general manner. 

We may simplify the determination of the Beta function parameters, 
avoiding the necessity to solve the cubic equation, (10), by the follow- 
ing procedure: 

Fix the variance of the symmetric Beta function, (m = .5), by 
ie tting 

d 2 (t* 

2 



p 



V 



P-0) 

~Tj 



(19) 



where d is an arbitrary positive constant — \jl2, or 



<f 2 ( x * 



m = 



5 



JL 

d2 



( 20 ) 



For the case where m = .5, equations (8), (9), and (20) yield the 
following: 
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( 21 ) 



CL = f r 



<i 2 - 12 



8 



Define C 




( 22 ) 



Then, for any value of m in its range, let ot and f be determined as 
follows : 



( = C , 0 — m — .5 

<X = C, .5 £ m l 



(23) 



The remaining parameter may be determined by equation (8). The variance 
of x*, as determined by equation (9), then varies with ra as shown in 



Figure 11. 




Change of Variance with Mode Position 



The mean, E [x*] , may be calculated directly by equation (12). The 
following close linear approximations for E[x*], for several values of d. 
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show the relative* weights given to (), M, and P in the determination of 
expected value, hence give an intuitive feeling for the relationship be- 
tween d and the amount of uncertainty with which m is located, 

d E |x*j (24) 

4 20 -f M 4 2P 

5 

6 0 4 4M + P 

6 

8 0 + 6M 4- P 

8 

By this relatively simple procedure, the parameters of the Beta 
function for any desired values of d may be determined for the purpose of 
the Monte Carlo calculations. The Probability distribution functions 
needed for the transformation from u* to x*, may be extracted from a table 
of the Incomplete Beta Function, fl] 
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APPENDIX III 



Use of the Triangular Probability Density in PERT 
Monte Carlo Calculations 

The triangular probability density function is naturally defined by 
the three time estimates, 0, M, and P. Since the Distribution function of 
this density nay be calculated by a simple formula, the transformation 
from a sample, u, drawn from a uniform (0,1) distribution, to a sample 
from the desired triangular density, may be accomplished by a single, 
simple calculation. When this procedure can be used, a considerable sav- 
ing in computer running time is effected. 

The technique is illustrated below for the random variable, y*> hav- 
ing the general triangular distribution with lower limit, 0, upper limit, 
P, and node, M, 

h 




Figure 15 

Triangular Density Function 



The density function of y* is: 






2(y-0) 
(P-0) (M-0) 



2(P- y) 

(P-0)(P-M) ’ 



0 , 



0 ^ y £ M 

M — y P 
elsewhere 



(25) 
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The distribution function, F^. v ( y ) z 



fy*(z)dz, is; 



v (y) 



(y-0)' 



(P-0)(M-0) 



0 ± y £ M 



( 26 ) 



- + M - 2PM + 2Py - y 2 , M ^ y ^ P 

P-0 (P-0) (p M) 

= 0 , y ^-0 

= 1, y^P 

Now lot u ; Py*(y)> and. solve for the inverse function, y=F v *^(u), 



and we got: 



y ; 0 + \/u(P-0)(H-0) , 0 £ u 6: 

P-0 



(27) 



= P - v(P-M) (P-0)(l-u) , ^ u ^1 

If we draw a random number, u, from the uniform (0,1) distribution, 
and determine y by means of the function, (27), then y* has the desired 
triangular distribution. 

The mean, E[y*3> is given exactly by: 

E[y*] = 0 + ^ f £ (28) 



2 

The variance, O' (y*), is given by: 



L *\ z 



ahy*) , 2^J . 



( 29 ) 



The terra on the right varies from 0, when M is at either extreme, to 



j P-0 

V6 



, when M is at the midpoint. 
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APPENDIX IV 



The SiiiCos Function 

Another Probability Function Suitable for PERT 
Monte Carlo Calculations 

Another probability function which is completely specified by its 
extremes and its mode is the following, which we shall call the SinCos 
function, for obvious reasons. 

The probability density function is: 



f 2 *(z> 



Tf 

2(P-0) 



sin 



7 T(z-0) 

2(M-0) 



t) £ z - M 



— . XL cos ^ - fe-M) . 
2(P -0) 2{P-M) 



M — z — P 



(30) 



a 0 , 



elsewhere 



The function has the following shape, where M can take any position be- 
tween and P. 







Figure 1G 

The SinCos Density Function 
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The distribution function, F ?v .(^) 



f *(t)dt, is 



J 



-Or 



p 2 *< z > - o. 

. Ml£ 

P-0 



l-cos JLSa-A). 
2(M-0) 



z ^ 0 



0 *£ z M 



( 31 ) 



!1'.Q 4 - P _M sin I r '.(.2.'il) , m ^ z -i P 
P-0 P-0 2(P-M) 



= 1 



z P 



Setting u -- p 2 * (z), and solving for z, we obtain: 



z : 0 -p 2(M-0) cos'^ 
TT 



1 - (P~0) u ] , 0 ^ u ^ M-0 

(M-0) | P-0 



= M 



, 2 < p -»> sin -i [t±L. itg’i 
+ rr L p ' M v ? ' 6 l 



(32) 



M-0 

P-0 



£ u 4. 1 



The above function transforms a random sample, u, from the uniform (0,1) 
distribution to the corresponding sample, z, from the SinCos distribution. 
The mean, E [ 2 !*], of the SinCos distribution is: 

(33) 



Elzi = ( TT-2)0 + (4-TT )M ± (T T -2 )P 

IT 



Considered as a weighted average, we sec that the extremes, 0 and P, re- 
ceive more weight than M. 

The variance, is; 



(j 2 (z *j . A( 1T - 3) (P-0 ) 2 + ( TT 2 -16T7V4O)(P-M) (M-0) 

In a form for numerical comparison with the PERT Beta variance: 
<f 2 (y.*) - 2.06(kl^ 2 - . 364 



(34) 



(35) 
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We see that 0 ^ ( z* ) is 1,88 i . 18K times the PERT Beta variance, P"^ ^ . 

6 

where - 1 K ^ 1, depending on the relative position of M. 

Although, with the SinCos function, the transformation from u* to 
z* can be accomplished by a single transformation equation, computer run- 
ning time may be as great, or greater, than the table look-up routine, 
due to the number of operations required to calculate the sin"** or the 
cos' 1 . 
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APPENDIX V 



FLOW CHART FOR RESOURCE ALLOCATION USING CRITICALITY 

INDEX 

• i 

For each T-C plan, form the following computer "word". 



| Activity No. 


ETB 


AT 


4C 
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